Crossover from Attractive to Repulsive Casimir Forces and Vice Versa 
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Systems described by an O(n) symmetrical cj> 4 Hamiltonian are considered in a d-dimensional 
film geometry at their bulk critical points. The critical Casimir forces between the film's boundary 
planes 23j, j = 1, 2, are investigated as functions of film thickness L for generic symmetry-preserving 
boundary conditions d n (f> ~ Cj(b. The L-dependent part of the reduced excess free energy per cross- 
sectional area takes the scaling form / rcs » D(ciL ,CiL '")/L d ~ l when d < 4, where a are 
scaling fields associated with the variables Cj, and $ is a surface crossover exponent. Explicit two- 
loop renormalization group results for the function D(ci,C2) at d — 4 — e dimensions are presented. 
These show that (i) the Casimir force can have either sign, depending on ci and C2, and (ii) for 
appropriate choices of the enhancements c,-, crossovers from attraction to repulsion and vice versa 
occur as L increases. 
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Macroscopic bodies that are immersed in a medium 
frequently experience long-range effective forces originat- 
ing from fluctuations in the medium. Such fluctuation- 
induced forces are ubiquitous in nature. A well-known 
example is the Casimir force between two metallic con- 
ducting plates caused by fluctuations of the electromag- 
netic field [l[ . Other important examples are the Casimir 
forces caused by confined thermal fluctuations, cither 
at critical points @, 0| or due to Goldstone modes [3]. 
Although predicted decades ago in a seminal paper by 
Fisher and de Gennes Q , such so-called "thermodynamic 
Casimir forces" were verified experimentally in a clear 
manner only recently — at first, indirectly by the thin- 
ning of He wetting layers near the lambda transition 
[f| and subsequently by their direct observation in bi- 
nary fluid mixtures Q. Current Monte Carlo simula- 
tions @, 0, H| were able to produce data in conformity 
with these experiments. 

Understanding fluctuation-induced forces is of great in- 
terest, for both technical and fundamental reasons, ft 
has recently been realized that quantum electromagnetic 
Casimir forces must be taken into account when design- 
ing micromechanical devices Analogous as well as 
thermodynamic Casimir forces are likely to be impor- 
tant also for microfluidic systems @, E3] ■ From a general 
vantage point, one of the most interesting aspects of fluc- 
tuation induced forces is their universality: They usually 
depend only on gross features of the medium, the macro- 
scopic bodies, and the geometry but are independent of 
microscopic details. 

A much studied case is the Casimir force Tq between 
two macroscopic parallel plates at a distance L, acting as 
boundaries of the medium in the z direction normal to the 
plates. It is frequently stated that, for a given medium 
and bulk dimension d, this force — and hence its sign — 
depends on the boundary conditions p on both boundary 
planes. Both the QED Casimir force and the thermody- 
namic Casimir force at a d-dimensional bulk critical point 



decay as inverse powers of L, the former (in three di- 
mension) as L~ 4 , the latter as L~ d . Their strengths are 
commonly characterized by dimensionless Casimir am- 
plitudes I\q , which are believed to be universal, though 
boundary condition-dependent 0, EH, E2, EH • For exam- 
ple, the critical Casimir force (measured in temperature 
units ksT and per cross-sectional area A) is convention- 
ally written as 



T c = -(d/dL)A^ ] L-^-^ = (d - 1) A^ p) L 



(1) 



If the boundary conditions are symmetric so that re- 
flection positivity holds, Tc is guaranteed to be attrac- 

< 0); for nonsymmetric 



14j (corresponding to A^ 



tive 

boundary conditions, repulsive Casimir forces may occur 
even for this simple slab geometry. 

The aim of this Letter is to show that the above pic- 
ture is oversimplified. Boundary conditions are scale- 
dependent properties. This entails that, even on length 
scales that are large compared to microscopic distances, 
the strengths of the critical Casimir forces cannot, in gen- 
eral, be characterized by constant universal amplitudes 
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i, . Rather, the get replaced by effective scale- 

(i.e., L-) dependent amplitudes. As L increases, they 
can change considerably. Particularly interesting is that 
even their signs may change so that originally attractive 
Casimir forces may turn repulsive as L increases and vice 
versa. Focusing on the case of critical Casimir forces, we 
shall present results for the associated scale-dependent 
amplitudes which show that — under appropriate con- 
ditions — smooth crossovers from repulsive to attractive 
as well as from attractive to repulsive Casimir forces are 
possible. Such crossovers should be accessible to experi- 
mental tests. 

To become more specific, let us consider an n- 
component <fi 4 theory on a slab 53 = M d_1 x [0, L] bounded 
in the z direction by a pair of planes Si at z = and QS2 
and z = L. For simplicity, we assume that these planes 
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do not give rise to interactions breaking the 0(n) sym- 
metry of the Hamiltonian TL. An appropriate choice then 
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where and are volume and surface integrals, re- 
spectively. This Hamiltonian is well known from the 
study of surface critical behavior. Let us recall some 
well-known facts needed below 

In Landau theory the Robin boundary conditions 



d n <p = Cj<fi on 05, 



(3) 



result, where d n means a derivative along the inner nor- 
mal. For 6j = oo they reduce to Dirichlct and for 6j — 
to Neumann boundary conditions. The physical signifi- 
cance of the interaction constants Cj is to account for local 
changes of the pair interactions near the planes 05 j . The 
larger the Cj, the stronger is the order parameter (p sup- 
pressed at 05.,. In Landau theory, Cj = corresponds to 
the special value Cj = c sp at which the plane 05., becomes 
critical exactly at the bulk transition temperature T c j,. 
More precisely, when Cj < c sp , a continuous phase transi- 
tion from a disordered phase to a (bulk-disordered) phase 
with long-range surface order at 05., occurs in the semi- 
infinite (L = oo) system at a temperature T C)S (cj) > T cb . 
The special value c sp of Cj at which T c>s (c sp ) = T c ^ spec- 
ifies a surface multicritical point, called "special" on the 
bulk critical line f = f c .fc. In Landau theory, both f C) b 
and c sp vanish. 

Beyond Landau theory, several important changes oc- 
cur. First, the boundary conditions © fluctuate — they 
hold in an operator sense, i.e., inside of averages [15J. 
Second, provided the multicritical point exists, — i.e., 
when d is sufficiently large that long-range surface order 
is possible at T > T c ^ — it gets shifted to nonzero val- 
ues of f Ct b and c sp , which depend on microscopical details 
(lattice constant a etc.). 

Thus, for critical enhancement Cj — c sp , a Robin 
boundary condition ([3]) with a nonuniversal c sp rather 
than a Neumann boundary condition applies (on the 
mesoscopic scale on which a continuum description is 
appropriate). This does not automatically rule out the 
validity of a Neumann boundary condition in the large- 
scale limit z — > oo with a <C z < £ (where £ is the bulk 
correlation length). The behavior of the order parame- 
ter near 05., follows from the boundary operator expan- 
sion (f>(x) ~ C(Az) (frlftj , where ^l^. is located on 05., 
at a distance Az from x. The short-distance behavior 
C(Az) ~ |Az|^! ~P>/ U is governed by the difference of 
the scaling dimensions (3/v and /3j p jv of <fi and 0|<8j> 
respectively. Only when their difference vanishes does a 
Neumann boundary condition hold on large length scales. 
While this is the case when d exceeds the upper critical 



dimension d* = 4 (since (3 = /3j p = 1/2 in Landau the- 
ory), it fails when d < d* because (3 > /3^ p . Thus, nei- 
ther on mesoscopic nor on large scales does a Neumann 
boundary condition hold at the special transition when 
d < d*. 

Conversely, one may choose Cj = 0, so that a Neu- 
mann boundary condition holds on a mesoscopic scale, 
and inquire again into the large-scale boundary condi- 
tions. Now Cj = translates into nonzero deviations 
from the multicritical point. Hence the 
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scaling fields Cj ~ 86j must vary under changes [i — > \il 
of the momentum scale. They become scale-dependent 
quantities Cj{€) that behave ~ l~^l v Cj near Cj = 0, 
where $ and v are the familiar surface crossover and 
bulk correlation length exponents, respectively. Depend- 
ing on whether Cj > or Cj < 0, they approach the 
fixed-point values c* rd = oo and c* x = — oo at which 
the fixed points describing the ordinary and extraordi- 
nary transitions of semi-infinite systems are located. The 
short-distance behaviors of <p near 25j are known in both 
cases. At the ordinary fixed point (c.j = +oo), one has 
cj) ~ z^-^^dn^l^, where /3f rd > (3 (for d < 4 and 
in Landau, theory where f3° rd = 1); at the extraordinary 
fixed point [c.j = -co), one has <fi ~ |Az|~ /3//,y . Hence, 
whenever the initial Cj > 0, Dirichlct boundary condi- 
tions hold asymptotically on large length scales at 05 j. 
The upshot is that, for generic values of Cj 6 (— oo, oo) 
and mesoscale boundary conditions ([3]), the boundary 
conditions of the full interacting theory will change under 
scale transformations, even in the Neumann case Cj = 0. 

To elucidate the consequences for the critical Casimir 
force, recall that the reduced free energy of the slab per 
cross-sectional area A — > oo can be decomposed as 



F/k B TA = Lf b + f s + f res (L) 



(4) 



into contributions from the bulk density the sur- 
face excess density f s , and an L-dependent residual part 
/ r es(i)- The behavior of these quantities at T c ^ can be 
analyzed via field-theoretic renormalization group (RG) 
methods. The RG equations satisfied by /b and f s upon 
renormalization at d < d* are inhomogeneous. How- 
ever, the one of f Tes (L) is known to be homogeneous 
U, 12, 3, EH- Solving it at T Ci t yields the scaling form 



/ ros (L)/n » D{ Cl L*l^c 2 L*l v ) 



(5) 



where the Cj now denote renormalized quantities Cj — 
H~ l Z^dtj involving a familiar renormalization factor Z c 
of Ref. [15(. The function D(ci,C2) is universal (up to 
nonuniversal metric factors). It replaces the amplitude 
A^/n in the first form of Eq. (fTJ), while the critical 



Casimir force becomes 



F c /nnV(c 1 L*/»,c 2 L*/ v )L- 



(6) 



with 



P(ci,c 2 )= [d-l + ($/ l /)(c 1 a ci +c 2 9 C2 )]^(c 1 ,c 2 ). (7) 
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We have computed the functions D and T> in d = 4 — e 
dimensions for general nonnegative values of Ci and C2 
to two-loop order, using e as a small parameter. The 
required free-energy terms involve summations over the 
spectrum {k m } of the operator — <9 2 on [0,L]. The dis- 
crete values k m are fixed by the boundary conditions ^ 
and depend on ci, c 2 , and L. We evaluated such mode 
sums by means of complex integration, employing a vari- 
ant of Abel-Plana techniques that facilitated the separa- 
tion of bulk and surface terms [l6| • 

To present our results, we introduce the functions 



Sd,c 2 (£) = In 



1 _ (ci -t)(c 2 -t) e _ 2t 

(Ci + t)(C2 + t) 



(8) 



j W r j-iyt^dt 

C1 ' C2 Jo (* + Cl) 2 (t + C 2 ) 2 e*-(t 2 -C?)(t»-C§)' 

(10) 

and the polynomials 



p(0,0) 

Ci,C 2 
Cl,C 2 

p(2,2) 

Cl,C 2 

p(l,0) 
ci,c 2 

p(2,l) 
ci,c 2 

(2,0) 
ci,c 2 



P 



2C{C%(C 1 + C 2 + C!C 2 ), 

2(c? + c3 + 2c 2 c 2 + 2 Cl c 2 + (c 2 + c 2 ) 2 ), 
2, 

2 Cl c 2 (c 2 + c 2 )( Cl + c 2 + cica) = Pj^. 

2(c 1 +c 2 + c 2 +c2)=P c ( i 1 ; c 2 2 ), 

2c lC2 (ci + c 2 + Cl c 2 ) = pC°£). (11) 



1 f 00 

A)(ci,c 2 ) = ^ J dtt 2 g c 



,(*), 



(9) Then our result for Z) can be written as 



D(ci,c 2 ) = D (ci,ca)+e^ 1- 



7 — ln7r 



1 f 00 

A)(ci,c 2 ) - y dtg CltC2 {t)t 2 \nt 



n + 2 
n + 8 



2 2 



<r,A=0 



+ o(e), (12) 



where 7 is the Euler-Mascheroni constant. 

To check this result, one can set (ci,c 2 ) to (00,00), 
(oo,0), and (0,0) and confirm by analytic calculation of 
the required integrals that the respective series (fT2| re- 
duce to the 0(e) results of Ref. 11] for the amplitudes 

A (ord,ord) /n) A (°^P) /nj and A (*P.sP)/ n . Note that the 

case (ci,c 2 ) = (0,0) is special: Unlike A^ rd ' ord ) and 
^(ord,sp)^ ^ e amplitude D(0,0) does not have an expan- 
sion in integer powers of e but involves also half-integer 
powers e fe / 2 , with k > 3 (besides powers of lne when 
fc>3)0. 

In Fig. [TJ we show a plot of Z?(ci,c 2 ) for the d = 3 
Ising case n = 1. It was obtained by numerical evalua- 
tion of the O(e) result (| 1 2[) at e = n = 1. As one sees, 
D changes sign along certain paths. The same is true 
for the scaling function ([7|), and hence for the critical 
Casimir force Tc- Moreover, such sign changes of Tc 
occur upon increasing L provided (ci, c 2 ) have appropri- 
ate values. That crossovers from attractive to repulsive 
Casimir forces and vice versa can occur is illustrated in 
Fig.H 

To put these results in perspective, consider the case 
b\ = 0. Here a Neumann boundary condition holds at 
2$i for the regularized theory on the mesoscopic scale on 
which the continuum description applies. One can de- 
rive this model from a simple cubic lattice spin model 
whose ferromagnetic nearest-neighbor bonds J xx i have 




FIG. 1: Scaling function D(ci,C2) for n = 1 and d = 3. 
To cover the full domain (0,oo) 2 , we plotted D(ci,C2) as a 
function of Cj/(1 +Cj). The zeros of D are depicted as thick 
lines. 

strengths Jj (J — 1,2) and J, depending on whether 
both sites x and x' belong to 03 j or at least one of them 
is not a boundary site. From the known (approximate) 
relation between Cj and Jj/ J [Hj], one sees that at d = 3 
the value {J\/J)q = 5/4 corresponds to cj = 0. This is 
less than the value p sp = 1.500(4) at which the multicrit- 
ical point of the d = 3 Ising model is located according to 
Monte Carlo simulations (l7|. Thus, a mesoscopic Neu- 
mann boundary condition at 5Bj corresponds to a sub- 
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FIG. 2: Scaled Casimir forces J~cL d as functions of L for 
(ci,c 2 ) = (0,0.1) (solid line) and (ci,c 2 ) = (10,0.1) (dashed 
line). In the first (second) case a crossover from attractive 
to repulsive (repulsive to attractive) Casimir forces occurs as 
L — > oo. 

critical enhancement cj > 0. As explained, a Dirichlet 
boundary condition applies in such a situation on large 
length scales. 

A similar crossover C2 = c\V^l v — > oo occurs also for 
the choice ox = 0.1 made in Fig. [2] In regimes where Ci 
and C2 are both small or both large, Tc must be attrac- 
tive. Yet in regimes where ci is sufficiently small while 
C2 is large, Tc is repulsive. Depending on our choices 
c\ = and c\ = 10 (3> C2), crossovers from attractive to 
repulsive Casimir forces and vice versa occur. 

These predictions should be testable by Monte Carlo 
simulations for three-dimensional Ising models of the 
kind described above and studied in Ref. [l7j |. Ideal 
experimental systems to measure the calculated scaling 
function would satisfy three criteria: (i) order-parameter 
dimension n = 1; (ii) non-symmetry-breaking bound- 
aries; (iii) tunability of the effective boundary pair in- 
teractions (ci and C2). In the case of 4 He at the lambda 
transition, the boundaries do not break the 0(2) sym- 
metry. However, a two-component order parameter is 
involved, a long-range ordered surface phase should not 
be possible at d — 3, and it is unclear to us how the pa- 
rameters 6j can be varied. Experimental studies of binary 
liquid mixtures seem to us a more promising alternative. 
For them, (i) is evidently satisfied. Since walls usually 
favor one or the other component, the Z2 symmetry is 
broken by linear boundary terms — Jig. hj(f>, where each 
hj can have either sign. It was demonstrated in Ref. [l8| 
that the values of these fields hj can be changed by chem- 
ically modifying the surface. It is also known that differ- 
ent signs of hj can be realized by proper choices of the 
mixtures and substrates [l£j. Hence it should be possi- 
ble to realize experimental setups with hi and ft-2 small 
and of opposite signs, or with hi m and /12 large and 
of equal signs. In both cases, sign-changing crossovers of 
the critical Casimir forces may be expected as L grows. 
The associated scaling functions would need separate cal- 



culations. 
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